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In the paper the closed Friedmann-Robertson-Walker model with quantization in the presence 
of the positive cosmological constant and radiation is studied. For analysis of tunneling probability 
for birth of an asymptotically deSitter, inflationary Universe as a function of the radiation energy 
a new definition of a "free" wave propagating inside strong fields is proposed. On such a basis, 
tunneling boundary condition is corrected, penetrability and reflection concerning to the barrier are 
calculated in fully quantum stationary approach. For the first time non-zero interference between 
the incident and reflected waves has been taken into account which turns out to play important 
role inside cosmological potentials and could be explained by non-locality of barriers in quantum 
mechanics. Inside whole region of energy of radiation the tunneling probability for the birth of 
the inflationary Universe is found to be close to its value obtained in semiclassical approach. The 
reflection from the barrier is determined for the first time (which is essentially differs on 1 at the 
£^ , energy of radiation close to the barrier height). The proposed method could be easily generalized 

on the cosmological models with the barriers of arbitrary shape, that has been demonstrated for 
the FRW-model with included Chaplygin gas. Result is stable for variations of the studied barriers, 
Q\ | accuracy are found to be 11-18 digits for all coefficients and energies below the barrier height. 
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W)| I. INTRODUCTION 

In order to understand what really happens in the formation of the Universe, many people came to the point 
of view that a quantum consideration of this process is necessary. After the publication of the first paper on the 
quantum description of Universe formation [14, 49], a lot of other papers appeared in this topic (for example, see 
Refs. [4, 19, 22, 37, 43-45, 52] and some discussions in Refs. [38, 46] with references therein). 

Today, among all variety of models one can select two approaches which are the prevailing ones: these are the 
• . Feynman formalism of path integrals in multidimensional spacetime, developed by the Cambridge group and other 
researchers, called the "Hartle-Hawking method" (for example, see Ref. [19]), and a method based on direct considera- 
£^ tion of tunneling in 4-dimensional Euclidian spacetime, called the "Vilenkin method" (for example, see Refs. [43-46]). 
■ Here, according to Ref. [47], in the quantum approach we have the following picture of the Universe creation: a closed 
Universe with a small size is formed from "nothing" (vacuum), where by the word "nothing" one refers to a quantum 
state without classical space and time. A wave function is used for a probabilistic description of the creation of the 
Universe and such a process is connected with transition of a wave through an effective barrier. Determination of 
penetrability of this barrier is a key point in the estimation of duration of the formation of the Universe, and the 
dynamics of its expansion in the first stage. 

However, in the majority of these models, with the exception of some exactly solvable models, tunneling is mainly 
studied in details in the semiclassical approximation (see Refs. [38, 46]). An attractive side of such an approach is 
its simplicity in the construction of decreasing and increasing partial solutions for the wave function in the tunneling 
region, the outgoing wave function in the external region, and the possibility to define and to estimate in an enough 
simply way the penetrability of the barrier, which can be used for obtaining the duration of the nucleation of the 
Universe. The tunneling boundary condition [46, 47] could seem to be the most natural and clear description, where 
the wave function should represent an outgoing wave only in the enough large value of the scale factor a. However, 
is really such a wave free in the asymptotic region? In order to draw attention on the increase of the modulus of 
the potential with increasing value of the scale factor a and increasing magnitude of the gradient of such a potential, 
acting on this wave "through the barrier" , then one come to a serious contradiction: the influence of the potential on 
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this wave increases strongly with a! Now a new question has appeared: what should the wave represent in general 
in the cosmological problem? This problem connects with another and more general one in quantum physics — the 
real importance to define a "free" wave inside strong fields. To this aim we need a mathematical stable tool to study 
it. It is unclear whether a connection between exact solutions for the wave function at turning point and "free" wave 
defined in the asymptotic region is correct. 

Note that the semiclassical formula of the penetrability of the barrier is constructed on the basis of 
wave which is defined concerning zero potential at infinity, i.e. this wave should be free outgoing in the 
asymptotic region. But in the cosmological problem we have opposite case, when the force acting on 
the wave increases up to infinity in the asymptotic region. At the same time, deformations of the shape 
of the potential outside the barrier cannot change the penetrability calculated in the framework of 
the semiclassical approach (up to the second order). An answer to such problem can be found in non- locality 
of definition of the penetrability in quantum mechanics, which is reduced to minimum in the semiclassical approach 
(i. e. this is so called "error" of the cosmological semiclassical approach). 

The problem of the correct definition of the wave in cosmology is reinforced else more, if one wants to calculate the 
incident and reflected waves in the internal region. Even with the known exact solution for the wave function there is 
uncertainty in determination of these waves! But, namely, the standard definition of the coefficients of penetrability 
and reflection is based on them. In particular, we have not found papers where the coefficient of reflection is defined 
and estimated in this problem (which differs essentially from unity at the energy of radiation close to the height of 
the barrier and, therefore, such a characteristics could be interesting from a physical point of view). Note that the 
semiclassical approximation put serious limits to the possibility of its definition at all [21]. 

Thus, in order to estimate probability of the formation of the Universe as accurately as possible, we need a fully 
quantum definition of the wave. Note that the non-semiclassical penetrability of the barrier in the cosmological 
problems has not been studied in detail and, therefore, a development of fully quantum methods for its estimation is 
a perspective task. 

Researches in this direction exist [2], and in these papers was estimated the penetrability on the basis of tunneling 
of wave packet through the barrier. However, a stationary boundary condition has uncertainty that could lead to 
different results in calculations of the penetrability. The stationary approach could allow to clarify this issue. It 
is able to give stable solutions for the wave function (and results in Ref. [29] have confirmed this at zero energy of 
radiation), using the standard definition of the coefficients of the penetrability and reflection, is more accurate to their 
estimation. 

Aims of this Chapter are: (1) to define the wave in the quantum cosmological problem; (2) to construct the fully 
quantum (non-semiclassical) stationary method of determination of the coefficients of penetrability of the barriers 
and reflection from them on the basis of such a definition of the wave; (3) to estimate how much the semiclassical 
approach differs in the estimation of the penetrability from the fully quantum one. In order to achieve this goal, we 
need to construct tools for calculation of partial solutions of the wave function. In order to resolve the questions 
pointed out above, we shall restrict ourselves to a simple cosmological model, where the potential has a barrier and 
internal above-barrier region. 



II. COSMOLOGICAL MODEL IN THE FRIEDMANN-ROBERTSON-WALKER METRIC WITH 

RADIATION 

A. Dynamics of Universe in the Friedmann— Robertson— Walker metric 

Let us consider a simple model of the homogeneous and isotropic Universe in Friedmann-Robertson-Walker (FRW) 
metric (see Ref. [48], p. 438; also see Rcfs. [9, 23, 40, 41]): 

ds 2 = -dt 2 + a 2 (t) ■ (J^-j + r 2 (d6 2 + si^edcj) 2 )^ , h(r) = 1 - kr 2 , (1) 

where t and r, 6, <f) are time and space spherical coordinates, the signature of the metric is (— , +, +, +) as in Ref. [41] 
(see p. 4), a(t) is an unknown function of time and k is a constant, the value of which equals +1, or —1, with 
appropriate choice of units for r. Further, we shall use the following system of units: U = c = 1. For k = —1, the 
space is infinite (Universe of open type), and for k = +1 the space is finite (Universe of closed type). For k = 1 one 
can describe the space as a sphere with radius a(t) embedded in a 4-dimensional Euclidian space. The function a(t) 
is referred to as the "radius of the Universe" and is called the cosmic scale factor. This function contains information 
of the dynamics of the expansion of the Universe, and therefore its determination is an actual task. 

One can find the function a(t) using the Einstein equations taking into account the cosmological constant A in this 
metric (we use the signs according to the chosen signature, as in Ref. [41] p. 8; the Greek symbols \i and v denote 
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any of the four coordinates t, r, 9 and </>): 

Rpv - \ 9^ R = 8-kG T„ v + A, (2) 

where i? M „ is the Ricci tensor, R is the scalar curvature, is the energy-momentum tensor, and G is Newton's 
constant. From (1) we find the Ricci tensor and the scalar curvature R: 



a 



ad d 2 h' 2a 2 + a'd + 2k 



Rtt — —3—, R rr — —r + 2 , . — . ., 

a h h hr 1 — kr A 

h'r { ' 

= R ef) sin 2 9, R 0e = adr 2 + 2d 2 r 2 - h - — + 1 = 2d 2 r 2 + ad 2 r 2 + 2kr 2 



D tt D , „rrp , oe D , 4>4> p 6 " + 6ad + 6fc 

i? = .9 R tt +9 R rr + 9 Ree+g R<t>4>= -j • (4) 

The energy-momentum tensor has the form (see [41], p. 8): T M „ = (p + p) U^U V + p g^ u , where p and p are 
energy density and pressure. Here, one needs to use the normalized vector of 4- velocity U* = 1, U r = U e = = 0. 
Substituting the previously calculated components (2) of the Ricci tensor i2„„, the scalar curvature (4), the components 
of the energy-momentum tensor T M „ and including the component p ra ,d(a), describing the radiation in the initial stage 
(equation of state for radiation: p(a) = /9 ra d(a)/3), into the Einstein's equation (2) at p = v = 0), we obtain the 
Friedmann equation with the cosmological constant (see p. 8 in Rcf. [41]; p. 3 in Ref. [9]; p. 2 in Rcf. [47]): 

^-^{Jt + ^H' «-5^, (5) 

where d is derivative a at time coordinate. From here, we write a general expression for the energy density: 



B. Action, lagrangian and quantization 

We define the action as in Ref. ([47], see (1), p. 2): 



U p\dx\ (7) 
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Substituting the scalar curvature (4), then integrating item at d by parts with respect to variable t, we obtain the 
lagrangian (see Ref. [47], (11), p. 4): 

£(a>)= 8^(^ 2 + fc -^ a2/3(a) )- (8) 

Considering the variables a and d as generalized coordinate and velocity respectively, we find a generalized momentum 
conjugate to a: 

d£(a,a) 3 . 

Pa = ^d— = -^G aa (9) 

and then hamiltonian: 

• n, -N 1 f 2lrG 2 2 3fc 4 1 

h(a,p ) = Pa-£(a,a) = -- j^-p a + a ^-^ - a p(a)j>. ( 10 ) 

The passage to the quantum description of the evolution of the Universe is obtained by the standard procedure of 
canonical quantization in the Dirac formalism for systems with constraints. In result, we obtain the Wheeler-De Witt 
(WDW) equation (see Ref. [47], (16)-(17), in p. 4, [14, 39, 49]), which can be written as 

{-£ + "<•>}*•>-*-*•>. "(«)=(ifc) 2t « 2 -^« 4 - *— 
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where ip(a) is wave function of Universe. This equation looks similar to the one-dimensional stationary Schrodinger 
equation on a semiaxis (of the variable a) at energy _E ra( j with potential V (a) . It is convenient to use system of units 
where 8nG = M~ 2 = 1, and to rewrite V (a) in a generalized form as 

V(a) = Aa 2 -Bo 4 . (12) 

In particular, for the Universe of the closed type (k = 1) we obtain A — 36, B = 12 A (this potential coincides with 
Ref. [2]). 



C. Potential close to the turning points: non-zero energy case 

In order to find the wave function we need to know the shape of the potential close to the turning points. Let us 
find the turning points at p ,i n and a tp , ut concerning the potential (12) at energy E^: 



a t P ,in = y^-y i-y 1 -—&-> at p-° ut = v2B"V V — ^ 

Let us expand the potential V(a) (13) in powers of g ou t = a — a tp (where the point at p .i n or atp.out is used as <2t p . 
Expansion is calculated at these points), where (for small q) we restrict ourselves to the liner term: 

V(q) = V Q + Vi q, (14) 

where the coefficients Vq and V\ are: 

V = V(a = o tp ,in) = V(a = a tPiOUt ) = Aa 2 p - B a 4 p = E rad , 



2B \ A 2 J \ V A 2 J 1 (15) 



Now eq. (15) transforms into a new form at variable q with potential V(q): 



dq 2 



-—rtq) + V iq v(q)=Q. (16) 



III. TUNNELING BOUNDARY CONDITION IN COSMOLOGY 



A. A problem of definition of "free" wave in cosmology and correction of the boundary condition 

Which boundary condition should be used to obtain a wave function that describes how the wave function leaves 
the barrier accurately? A little variation of the boundary condition leads to change of the fluxes concerning the barrier 
and, as result, it changes the coefficients of penetrability and reflection. So, a proper choice of the boundary condition 
is extremely important. However before, let us analyze how much the choice of the boundary condition is natural in 
the asymptotic region. 

• In description of collisions and decay in nuclear and atomic physics potentials of interactions tend to zero 
asymptotically. So, in these calculations, the boundary conditions are imposed on the wave function 
at infinity. In cosmology we deal with another, different type of potential: its modulus increases with increasing 
of the scale factor a. The gradient of the potential also increases. Therefore, here there is nothing common 
to the free propagation of the wave in the asymptotic region. Thus, a direct passage of the application of the 
boundary condition in the asymptotic region into cosmological problems looks questionable. 

• The results in Ref. [29], which show that when the scale factor a increases the region containing solutions for 
the wave function enlarges (and its two partial solutions), reinforce the seriousness of this problem. According 
to Ref. [29], the scale factor a in the external region is larger, the period of oscillations of each partial solution for 
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the wave function is smaller . One has to decrease the time-step and as a consequence increase the calculation 
time. This increases errors in computer calculations of the wave function close the barrier (if 
it is previously fixed by the boundary condition in the asymptotic region). From here a natural 
conclusion follows on the impossibility to use practically the boundary condition at infinity for calculation of 
the wave (in supposition if we know it maximally accurately in the asymptotic region), if we like to pass from 
the semiclassical approach to the fully quantum one. Another case exists in problems of decay in nuclear and 
atomic physics where calculations of the wave in the asymptotic region are the most stable and accurate. 

• One can add a fact that it has not been known yet whether the Universe expands at extremely large value of the 
scale factor a. Just the contrary, it would like to clarify this from a solution of the problem, however imposing 
a condition that the Universe expands in the initial stage. 

On such a basis, we shall introduce the following definition of the boundary condition: 

The boundary condition should be imposed on the wave function at such value of the scale 
factor a, where the potential minimally acts on the wave, determined by this wave function. 

The propagation of the wave defined in such a way is close to free one for the potential and at used value of the scale 
factor a (we call such a wave conditionally "free"). However, when we want to give a mathematical formulation of 
this definition we have to answer two questions: 

1. What should the free wave represent in a field of a cosmological potential of arbitrary shape? How could it be 
defined in a correct way close to an arbitrary selected point? 

2. Where should the boundary condition be imposed? 

To start with, let us consider the second question namely where we must impose the boundary condition on the 
wave function. One can suppose that this coordinate could be (1) a turning point (where the potential 
coincides with energy of radiation), or (2) a point where a gradient from the potential (having a sense 
of force of interaction) becomes zero, or (3) a point where the potential becomes zero. But the clear 
condition of free propagation of the wave is the minimal influence of the potential on this wave. So, 
we define these coordinate and force in the following way: 

The point in which we impose the boundary condition is the coordinate where the force acting 
on the wave is minimal. The force is defined as minus the gradient of the potential. 

It turns out that according to such a (local) definition the force is minimal at the external turning point a tPiOUt . 
Also, the force, acting on the wave incident on the barrier in the internal region and on the wave reflected from it, has 
a minimum at the internal turning point a tPi i n . Thus, we have just obtain the internal and external turning points 
where we should impose the boundary conditions in order to determine the waves. 

B. Boundary condition at a = 0: stationary approach versus non-stationary one 

A choice of the proper boundary condition imposed on the wave function is directly connected with the question: 
could the wave function be defined at a = 0, and which value should it be equal to at this point in such a case? 
The wave function is constructed on the basis of its two partial solutions which should be linearly independent. In 
particular, these two partial solutions can be real (not complex), without any decrease of accuracy in determination 
of the total wave function. For any desirable boundary condition imposed on the total wave function, such methods 
should work. In order to achieve the maximal linear independence between two partial solutions, we choose one 
solution to be increasing in the region of tunneling and another one to be decreasing in this tunneling region. For the 
increasing partial solution we use as starting point a x the internal turning point a tp . i n at _E ra d ^ or zero point 
a x = at EVad = 0. For the second decreasing partial solution the starting point a x is chosen as the external turning 
point a tp , out- Such a choice of starting points turns out to give us higher accuracy in calculations of the total wave 
function than starting calculations of both partial solutions from zero or from only one turning point. 

In order to obtain the total wave function, we need to connect two partial solutions using one boundary condition, 
which should be obtained from physical motivations. According to analysis in Introduction and previous 
section, it is natural not to define the wave function at zero (or at infinity), but to find outgoing wave 
at finite value of a in the external region, where this wave corresponds to observed Universe at present 
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time. But, in practical calculations, we shall define such a wave at point where forces minimally act 
on it. This is an initial condition imposed on the outgoing wave in the external region 1 . 

Let us analyze a question: which value has the wave function at a = 0? In the paper the following ideas are used: 

• the wave function should be continuous in the whole spatial region of its definition, 

• we have outgoing non-zero flux in the asymptotic region defined on the basis of the total wave function, 

• we consider the case when this flux is constant. 

The non-zero outgoing flux defined at arbitrary point requires the wave function to be complex and non-zero. The 
condition of continuity of this flux in the whole region of definition of the wave function requires this wave function to 
be complex and non-zero in the entire region. If we include point a = into the studied region, then we should obtain 
the non-zero and complex wave function also at such point. If we use the above ideas, then we cannot obtain zero 
wave function at a = 0. One can use notions of nuclear physics, field in which the study of such questions and their 
possible solutions have longer history then in quantum cosmology. As example, one can consider elastic scattering of 
particles on nucleus (where we have zero radial wave function at r = 0, and we have no divergences), and alpha decay 
of nucleus (where we cannot obtain zero wave function at r = 0) . A possible divergence of the radial wave function 
at zero in nuclear decay problem could be explained by existence of source at a point which creates the outgoing flux 
in the asymptotic region (and is the source of this flux). Now the picture becomes clearer: any quantum decay could 
be connected with source at zero. This is why the vanishing of the total wave function at a = 0, after introduction of 
the wall at this point (like in Ref. [2]), is not obvious and is only one of the possibilities. 
If we wanted to study physics at zero a = 0, we should come to two cases: 

• If we include the zero point into the region of consideration, we shall use to quantum mechanics with included 
sources. In such a case, the condition of constant flux is broken. But a more general integral formula of non- 
stationary dependence of the fluxes on probability can include possible sources and put them into frameworks of 
the standard quantum mechanics also (sec cq. (19.5) in Ref. [21], p. 80). Perhaps, black hole could be another 
example of quantum mechanics with sources and sinks. 

• We can consider only quantum mechanics with constant fluxes and without sources. Then we should eliminate 
the zero point a = from the region of our consideration. In this way, the formalism proposed in this paper 
works and is able to calculate the penetrability and reflection coefficients without any lost of accuracy. 

This could be a stationary picture of interdependence between the wave function at zero and the outgoing flux in the 
asymptotic region. In order to study the non-stationary case, then we need initial conditions which should define also 
the evolution of the Universe. In such a case, after defining the initial state (through set of parameters) it is possible 
to connect zero value of wave packet at a = (i. e. without singularity at such a point) with non-zero outgoing flux 
in the asymptotic region. In such direction different proposals have been made in frameworks of semiclassical models 
in order to describe inflation, to introduce time or to analyze dynamics of studied quantum system (for example, see 
[17, 42]). 

IV. DIRECT FULLY QUANTUM METHOD 

A. Wave function of Universe: calculations and analysis 

The wave function is known to oscillate above the barrier and increase (or decrease) under the barrier without 
any oscillations. So, in order to provide an effective linear independence between two partial solutions for the wave 
function, we look for a first partial solution increasing in the region of tunneling and a second one decreasing in this 
tunneling region. To start with, we define each partial solution and its derivative at a selected starting point, and 
then we calculate them in the region close enough to this point using the method of beginning of the solution presented 
in Subsection IV D 1. Here, for the partial solution which increases in the barrier region, as starting point we use the 
internal turning point at p ,i n at non-zero energy i? ra d or equals to zero a — at null energy -E ra d, and for the second 
partial solution, which decreases in the barrier region, we choose the starting point to be equal to the external turning 
point a tp , out- Then we calculate both partial solutions and their derivatives in the whole required range of a using the 



For example, on the basis of such a boundary condition for a-decay problem we obtain the asymptotic region where the wave function 
is spherical outgoing wave. 
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method of continuation of the solution presented in Subsection IV D 2, which is improvement of the Numerov method 
with constant step. In this way, we obtain two partial solutions for the wave function and their derivatives in the 
whole studied region. 

In order to clarify how the proposed approach gives convergent (stable) solutions, we compare our results with the 
paper of [2]. Let us consider the behavior of the wave function. The first partial solution for the wave function 
and its derivative in my calculation are presented in Fig. 1, which increase in the tunneling region and have been 
obtained at different values of the energy of radiation -E ra d. From these figures one can see that the wave function 




Parameter a Parameter a Parameter a 

FIG. 1: The first partial solution for the wave function and its derivative at different values of the energy of radiation E TBj< x, 
increasing in the tunneling region. The blue plot represents the wave function; the green one, the derivative of this wave 
function): (a) £ rad = 10; (b) £ rad = 1000; (c) £ rad = 2000 

satisfies the rules satisfied by the wave function inside the sub-barrier and in above-barrier regions [51]. Starting from 
very small a, the wave function has oscillations and its maxima increase monotonously with increasing of a. This 
corresponds to the behavior of the wave function in the internal region before the barrier (this becomes more obvious 
after essential increasing of scale, see left panel in Fig. 2). Moreover, for larger values of a, the wave function increases 
monotonously without any oscillation, that points out on transition into the tunneling region (one can see this in a 
logarithmic presentation of the wave function, see central panel in Fig. 2). A boundary of such a transformation in 
behavior of the wave function must be the point of penetration of the wave into the barrier, i. e. the internal turning 
point a tp! i n . Further, with increasing of a the oscillations are appeared in the wave function, which could be possible 
inside the above barrier region only (in the right panel of Fig. 2 one can see that such a transition is extremely 
smooth that characterizes the accuracy of the method positively). The boundary of such a transformation in 
the behavior of the wave function should be the external turning point Otp.out- Like Ref. [29], but at 
arbitrary non-zero energy _E la d we obtain monotonous increasing of maximums of the derivative of the 
wave function and smooth decreasing of this wave function in the external region. One can see that 
the derivative is larger than the wave function. At large values of a we obtain the smooth continuous 
solutions up to a — 100 (in Ref. [2] the maximal presented limit is a = 30). 




Parameter a Parameter a Parameter a 



FIG. 2: The first partial solution for the wave function and its derivative at the energy of radiation £Wd = 2000. The blue line 
represents the wave function; the green one, the derivative of this wave function) 

In Fig. 3, it is presented the second partial solution of the wave function and its derivative at different values of the 
energy of radiation According to the analysis, this solution close to the turning points, in the tunneling region, 
in the sub-barrier and above-barrier regions looks like the first partial solution, but with the difference that now the 
maxima of the wave function and their derivatives are larger essentially in the external region in a comparison with 
the internal region, and amplitudes in the tunneling region decrease monotonously. 

Comparing the previous pictures of the wave function with the results of Ref. [2] , one can see that the wave function, 
in this approach, is essentially more continuous, has no divergencies and its behavior is everywhere clear. From here we 
conclude that the developed method for the determination of the wave function and its derivative at arbitrary energy of 
radiation is essentially more quick, more stable and accurate in comparison with the non- stationary quantum approach 
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Parameter a 



Parameter a 



Parameter a 



FIG. 3: The second partial solution for the wave function and its derivative at different values of the energy of radiation 
i?rad, decreasing in the tunneling region (curve 1, blue, is for the wave function; curve 2, green, for the derivative of this wave 
function): (a) E Tad = 10; (b) E rad = 1000; (c) £ rad = 2000 



in Ref. [2]. Note that: 

• With increasing a, the period of the oscillations, both for the wave function and its derivative, decreases uniformly 
in the external region and increases uniformly in the internal region (this result was partially obtained earlier 
in Ref. [29] at E Iad = 0). 

• At larger distance from the barrier (i. e. for increasing values of a, in the external region, and at decreasing 
value of a, in the internal region) it becomes more difficult to get the convergent continuous solutions for the 
wave function and its derivative (this result was partially obtained earlier in Ref. [29] at E la a = 0). 

• A number of oscillations of the wave function in the internal region increases with increasing of the energy of 
radiation E Tad (this is a new result). 



B. Definition of the wave minimally interacting with the potential 

Now we shall be looking for a form of the wave function in the external region, which describes accurately the wave, 
whose propagation is the closest to the "free" one in the external region at the turning point atp.out and is directed 
outside. Let us return back to eq. (16) where the variable q = a — a tPi0 ut has been introduced. Changing this variable 
to 

£= |^ (out) | 1/3 <7, (17) 

this equation is transformed into 

|^(£)+^(0=0. (18) 

From quantum mechanics we know two linearly independent exact solutions for the function in this equation 
— these are the Airy functions Ai(£) and Bi(£) (see Ref. [1], p. 264-272, 291-294). Expansions of these functions 
into power series at small £, their asymptotic expansions at large |£|, their representations through Bessel functions, 
zeroes and their asymptotic expansions are known. We have some integrals of these functions, and also the form of 
the Airy functions in the semiclassical approximation (which can be applied at large |£|). In some problems of the 
analysis of finite solutions in the whole range of £ it is convenient to use the integral representations of the Airy 
functions (see eq. (10.4.32) in Ref. [1], p. 265. In eq. (10.4.1) we took into account the sign and a = 1/3): 

+ OG 

Ai(±0 = I J 



Bi( ± f) . I / 

o 

Furthermore, we shall be interested in the solution <p(£) which describes the outgoing wave in the range of a close 
to the dtp point. However, it is not clear what the wave represents in general near the point at p , and which linear 
combination of the Ai (£) and Bi (£) functions defines it in the most accurate way. 



cos [ — =F £, u ) du, 



exp 



-y Ttu) +sin T^u 



(19) 



du. 
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The clearest and most natural understanding of the outgoing wave is given by the semiclassical consideration of the 
tunneling process. However, at the given potential the semiclassical approach allows us to define the outgoing wave 
in the asymptotic region only (while we can join solutions in the proximity of at p by the Airy functions). But it is not 
clear whether the wave, defined in the asymptotic region, remains outgoing near the a tp - During the whole path of 
propagation outside the barrier the wave interacts with the potential, and this must inevitably lead to a deformation 
of its shape (like to appearance of a phase shift in the scattering of a wave by a radial potential caused by interaction 
in scattering theory). Does the cosmological potentials deform the wave more than the potentials used 
for description of nuclear collisions in scattering theory? Moreover, for the given potential there is a problem 
in obtaining the convergence in the calculation of the partial solutions for the wave function in the asymptotic region. 
According to our calculations, a small change of the range of the definition of the wave in the asymptotic region leads 
to a significant increase of errors, which requires one to increase the accuracy of the calculations. Therefore, we shall 
be looking for a way of defining the outgoing wave not in the asymptotic region, but in the closest vicinity of the point 
of escape, a tp - In a search of solutions close to the point a tp , i. e. at small enough the validity of the semiclassical 
method breaks down as |£| approaches zero. Therefore, we shall not use the semiclassical approach in this paper. 

Assuming the potential V(a) to have an arbitrary form, we define the wave at the point a tp in the following way. 



Definition 1 (strict definition of the wave). The wave is a linear combination of two partial solutions of the wave 
function such that the change of the modulus p of this wave function is approximately constant under variation of a: 

-> 0. (20) 

a=a tp 

According to this definition, the real and imaginary parts of the total wave function have the closest behaviors 
under the same variation of a, and the difference between possible maximums and minimums of the modulus of the 
total wave function is the smallest. For some types of potentials (in particular, for a rectangular barrier) it is more 
convenient to define the wave less strongly. 




Definition 2 (weak definition of wave): 

The wave is a linear combination of two partial solutions of wave function such that the modulus p changes 
minimally under variation of a: 



Ta P(a) 



-> 0. (21) 

a—a tp 



According to this definition, the change of the wave function caused by variation of a is characterized mainly by its 
phase (which can characterize the interaction between the wave and the potential). 
Subject to this requirement, we shall look for a solution in the following form: 

p(0=T-* (+) (0, (22) 

where 

Umax 

P 3 

= J exp±i(-y + /(0u)d«. (23) 
o 

where T is an unknown normalization factor, /(£) is an unknown continuous function satisfying /(£) — > const at 
£ — > 0, and u ma x is the unknown upper limit of integration. In such a solution, the real part of the function /(£) gives 
a contribution to the phase of the integrand function, while the imaginary part of /(£) deforms its modulus. 

Let us find the first and second derivatives of the function (a prime denotes a derivative with respect to £): 



"max g 

^* (±) (0 = ±we J /'uexp±i(-y + /(£)«) du, 



Umax 

^* (±) (0 = / (± if"u-{f) V) exp±i(-^ + /(0«)d«. 



(24) 
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From this we obtain: 

72 



^* (±) (0+e* (±) (0= / (± */"«-(/') 2 « 2 +e) exp±i(-y+/(0t 



Considering the solutions at small enough values of |£|, we represent /(£) in the form of a power series: 



+ OC 



n=0 

where /„ are constant coefficients. The first and second derivatives of /(£) are 

, +00 +00 

a 



/'(o = ^ /(o - E n /» = E( n + !) 



n=l 

+ 00 



/"(o = ^2 /(o = E( n + !) ( n + 2 ) 

? n=0 

Substituting these solutions into eq. (24), we obtain 

«max ( 

!_ * (±) (0 +f *<±> (0 = J (± 2iu /2 - u * ^ + 

/ \ +°° r 

+ ± 6iu h - 4« 2 A/2 + 1) C + + E ± it* (n + l)(n + 2) / n+2 - 



-u 2 £ (n - m + l)(m + 1) /„_ m +i/ m +i 



m=0 



; n | exp±i(-y + /t 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



Considering this expression at small |£|, we use the following approximation: 

exp±i(-y + /u) ->■ exp±z(-y + f u). 
Then from eq. (18) we obtain the following condition for the unknown /„: 

J (± 2iu f 2 - u 2 / 2 ) cxp±z (- y + /ow) + 


W m ax g 

+ e • y (± 6i« / 3 - 4u 2 A A + l) cxp ±» (~ + /o«) du + 


+ oo U J\ JL 

+ Ef' / ±Mi(n+l)(n + 2)/„ +2 - M 2 ^(n-m+l)(m+l)/„_ m+1 / ra+1 
™= 2 m=0 

x cxp±i (-y + f a uj du = 0. 
Requiring that this condition is satisfied for different £ and with different powers n, we obtain the following system: 

C° : / (± h - u 2 A 2 ) exp ±i (- y + /„«) du = 0, 



1i max g 

: y (± 6zw / 3 - 4« 2 A/2 + l) cxp ±i (- y + Am) d« - 0, 



(31) 



± m (n + l)(n + 2) / n+2 - w 2 E ( n _ m + 1 )( m + x ) A-m+i/m+i 



m=0 
~3 



x cxp ±i ^ — — + A u ) rfu = 0. 
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Assuming that the coefficients f and f\ are known, we find the following solutions for the unknown / 2 , f-$ and /„: 

f 2 7 (±) M 4A f (±) / (±) - 7 (±) 
f (±) = ± LL . - 2 — = ± 4/1/2 J ° , (32) 

2 2z j(±)' 7 3 6zJ (±) 

£ (n - m + l)(m + 1) /^ m +i /i+i r(±) , ooN 

»(±) _ m=0 Jj. (33) 

/n+2 ~ i(n + l)(n + 2) "j(±)' 

where the following notations for the integrals have been introduced: 

Umax Umax 

f 3 r 3 

^o ±) = / exp±i (--y + /o^ du, = u exp±i (-y + fou) du, (34) 



,3 



Mmax / U J 

/ AH-y + ^/d, (35) 


Thus, we see that the solution (22) taking into account eq. (23) for the function y (£) has arbitrariness in the choice 
of the unknown coefficients fo, fx and the upper limit of integration, u max . However, the solutions found, eqs. (32), 
define the function /(£) so as to ensure that the equality (22) is exactly satisfied in the region of a close to the escape 
point dtp- This proves that the function tp(£) in the form (22), taking into account eq. (23) for an arbitrary choice 
of fo, fx and w max is the exact solution of the Schrodinger equation near the escape point a tp . In order to write the 
solution ^(£) in terms of the well-known Airy functions, Ai (£) and Bi (£), we choose 

fo = 0, fx = 1. (36) 

For such a choice of the coefficients fo and fx, the integrand function in the solution (23) (up to £ 2 ) has a constant 
modulus and a varying phase. Therefore, one can expect that the solution (22) at the turning point a tp describes the 
wave accurately. 



C. Total wave function 



Having obtained two linearly independent partial solutions <fi(a) and </?2(«), we can write the general solution (a 
prime is for the derivative with respect to a) as: 

<p(a) = T-(C 1 <p 1 (a) + C2<to(a)), (37) 

tf ^2 - *'<p2 
°1 — 7 7 

where T is a normalization factor, Cx and C2 are complex constants found from the boundary condition introduced 
above: the f>{a) function should represent an outgoing wave at turning point a tp , ut- Fig. 4 plots the total wave 
function calculated in this way for the potential (12) with parameters A = 36, B = 12 A at A = 0.01 at different 
values of the energy of radiation i? ra d . One can see that the number of oscillations of the wave function in the internal 
region increases with increasing of the energy of radiation. Another interesting property are the larger maxima of the 
wave function in the internal region at smaller distances to the barrier for arbitrary energy (result found for the first 
time). 

In Fig. 5 it has been shown how the modulus of this wave function changes at selected values of the energy of 
radiation. From these figures it becomes clear why the coefficient of penetrability of the barrier is extremely small 
(up to the energy E la( x = 2000). In order to estimate, how effective is the boundary condition introduced above in 
building up the wave on the basis of the total wave function close to the external turning point a tp , ut, it is useful to 
see how the modulus of this wave function changes close to this point. In Fig. 6 we plot the modulus of the found 
wave function close to the turning points at the energy of radiation E ra( x = 2000 is shown. Here, one can see that the 
modulus at a tp ,out is practically constant (see left panel in Fig. 6). It is interesting to note that the modulus of the 
wave function, previously defined, does not change close to the internal turning point atp, in, and is close to maximum 
(see right panel in Fig. 6). 



C 2 = 



(38) 



CI, O.l-r, ,-..< + 
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Parameter a 



FIG. 4: The wave function at selected values of the energy of radiation i5 ra d (the blue line, represents the real part of the wave 
function; the green line the imaginary part of the wave function): (a) E la _d = 10; (b) -E ra d = 1000; (c) E la _d = 2000 
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FIG. 5: The behavior of the modulus of the wave function at the selected energies of radiation E la ,d'- (a) U ra d = 10; (b) 
Brad = 1000; (c) £ rad = 2000. 



D. Calculations of the wave function of Universe 

1. Method of calculations of the wave function close to an arbitrary selected point a x 

Here, we look for the regular partial solution of the wave function close to an arbitrary selected point a x . Let us 
write the wave function in the form: 

ip(a) — c 2 K (a - a x ) n = c 2 J2 bn a™ , a — a - a x (39) 

n=0 n=0 

and rewrite the potential through the variable a: 

V(a) = C + C\a + C 2 a 2 + C 3 a 3 + Cia 4 , (40) 

where 

C = Aal-Bai, 

Ci = 2a x (A-Bal)-2Bal=2Aa x -4Bal, 

C 2 = A- Ba 2 x -ABa 2 x - B a 2 x = A-6Ba 2 x , (41) 
C 3 = -2Ba x - 2Ba x = -ABa x , 
C 4 = -B. 

Substituting the wave function (39), its second derivative and the potential (40) into Schrodinger equation, we obtain 
recurrent relations for unknown b n : 

, (Co-E)bo , (Co - E) bi + Cx b (C - E) b 2 + Cy h + C 2 b 

b 2 g ' 63 6 ' 64 = 12 ' [ ' 

, _ (Co - E) b 3 + Ci b 2 + C 2 b 1 + C 3 b m 
5 20 ' 

(Co — E) b n + C\ 6 n _i + C2 6 n _2 + C3 b n s + C4 6 n _4 

b n+2 = -, — — TT7 — -rr at n > 4. (44 

(n + 1) (n + 2) 
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FIG. 6: The behavior of the modulus of the total wave function at the energy of radiation -E ra d = 2000, close to the turning 
points (for a tp ,i n = 8.58, a tp , ut = 15.04, see also Table 1): (a) the modulus decreases monotonously in the tunneling region, 
with increasing of a. It shows maxima and holes connected with the oscillations of the wave function in the external region, but 
the modulus is not equal to zero (thispoints out the existence of a non-zero flux); (b) when a increases, the modulus reaches a 
minimum close to the external turning point dtp, out (this demonstrates the practical fulfillment of the definition for the wave 
at such a point); (c) transition close to a tp ,m is shown, where at increasing of a the modulus with maximums and holes is 
transformed rapidly into a monotonously decreasing function without maximums and holes. This is connected with transition 
to the region of tunneling. 



Given the values of bo and b\ and using eqs. (42)-(44) one can calculate all b n needed. At limit £7 ra( j — > and at 
a x — all found solutions for bi transform into the corresponding solutions (40), early obtained in [29] at -E ra d = 0. 
Using C2 = 1, from eqs. (39) we find: 

b = tp(a x ), h = (p'(a x ). (45) 

So, on the basis of the coefficients bo and b\ one can obtain the values of the wave function and its derivative at point 
a x . Imposing two different boundary conditions via bo and 6i, we obtain two linear independent partial 
solutions <pi(a) and ^2(0) f° r the wave function. Using the internal turning point at Pj i n as the starting point, 
we calculate the first partial solution which increases in the barrier region (we choose: bo = 0.1, b\ = 1), and using 
the external turning point dtp, out as the starting point, we calculate the second partial solution which decreases in 
the barrier region (we choose: 60 = 1> b\ = —0.1). Such a choice provides effectively a linear independence between 
two partial solutions. 



2. Method of continuation of the solution 
Let us rewrite equation (18) in such a form 2 : 

<p"(a) = f(a)<p(a). (46) 

Let {a ra } be a set of equidistant points a n = ao + nh. Denoting the values of the wave function ip (a) at points a n as 
ip n , we have constructed an algorithm of the ninth order to determine (p n +l an d ^>' n when ip n and ip n -i are known: 

901 - ffll 301 ~ 311 / 47 \ 

> 2 , 3io ™ 3oo , ^ (u o\ 

3oi - 3n 3oi - 3n 



2 Here, we used the algorithm of [50] 
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whore 

goo = 2 + ^/„ + |/ l 4 (/: + /2) + |^(/( 4 )+4(/;) 2 + 7/„.C + /„ 3 ) + 

+ I h& + 16 fn ; " 4) + 26 & ; " 3) + 15 ( / «) 2 + 22 ^ ^ + 28 fn ( / ") 2 + / «) ' 

.901 = 1^2/; + !^ (4 + 6 /„ /;) + | h* (6 /( 5 ) + 24 /„ /( 3 ) + 48 f'Jl + 12 /* , (48) 



5io = | h 3 f' n + | h 5 (f& + 4 /„ /;) + | h 7 (/W + 11 /„ /( 3 ) + 15 /X + 9 /* /;) , 
9ii = 2/ l +|/, 3 /ri + |/ l 5 (3/: + /2) + |/ l 7 (5/( 4 ) + 13/„/: + 10(/;) 2 + /„ 3 ). 
A local error of these formulas at point a n equals to: 

5n = ^h 10 f^\ (49) 

E. The penetrability and reflection in the fully quantum approach 

Let us analyze whether a known wave function in the whole region of its definition allows us to determine uniquely 
the coefficients of penetrability and reflection. 

1. Problem of interference between the incident and reflected waves 

Rewriting the wave function </? t otai in the internal region through a summation of incident </?; nc wave and reflected 
<p Te { wave: 

'/'total = </>inc + tp le i, (50) 

we consider the total flux: 



j ('/'total) 

where 



(<Anc + ¥>ref)v(<pf nc + ^ ef ) - h. C.) 



(51) 



jinc = i (vincVv?* nc - h. C.) , j rc f = 1 (</?ref V</J* ef - h. C.) , 

jmixed = i (<PincV<Pre f + if Ie { V<p* nc - h. C.) . 



(52) 



The jmixod component describes interference between the incident and reflected waves in the internal region (let us 
call it mixed component of the total flux or simply flux of mixing). From the constancy of the total flux jtotai we find 
the flux j tr for the wave transmitted through the barrier, and: 

jinc = jtr - jrcf - jmixod, jtr = jtotai = COnSt. (53) 

Now one can see that the mixed flux introduces ambiguity in the determination of the penetrability and reflection for 
the same known wave function. 



F. Determination of the penetrability, reflection and interference coefficients 

In quantum mechanics the coefficients of penetrability and reflection are defined considering the potential as a whole, 
including asymptotic regions. However, in the radial calculation of quantum decay such a consideration 
depends on how the incident and reflected waves are defined inside finite internal region from the left 
of the barrier. The question is: does the location of such a region influence the penetrability and 
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reflection? In order to obtain these coefficients, we shall include into definitions coordinates where the fluxes are 
defined (denote them as x\ e f t and alight ) : 

rp f n _ jtr (bright ) „, ■> _ Jrcf (^lcft) , ,/ s _ jmixed (^lcft ) (t . ,n 

-L (Xleft, bright) — - — 7 r, n{Xi ett ) — -, M(X\ c{t )-— — — . (54) 

Jinc (^left ) Jinc (^left ) Jinc (#left j 

So, the T and R coefficients determine the probability of transmission (or tunneling) and reflection of the wave rela- 
tively the region of the potential with arbitrary selected boundaries xi f t , x r ; g ht- When x r i g ht tends to the asymptotic 
limit, the coefficient defined before should transform into standard ones. From eqs. (53) and (54) we obtain (j tr and 
jref are directed in opposite directions, j- lnc and jt r — in the same directions): 

\T\ + \R\ - M = 1. (55) 

Now we see that the condition \T\ + \R\ = 1 has sense in quantum mechanics only if there is no interference between 
incident and reflected waves, and for this is enough that: 

jmixed = 0. (56) 

A new question appears: does this condition allow to separate the total wave function into the incident and reflected 
components in a unique way? It turns out that the choice of the incident and reflected waves has essential influence on 
the barrier penetrability, and different forms of the incident </?i nc and reflected <p re f waves can give zero flux j m j x . Going 
from the rectangular internal well to the fully quantum treatment of the problem would become more complicated. 



G. Wave incident on the barrier and wave reflected from it in the internal region 

One can define the incident wave to be proportional to the function , 3>( + ) and the reflected wave to be proportional 
to the function 

(/'total (a) = <Anc ( a ) + ¥>ref (a), fine (a) = we • *(+) (a), Lp rcf (a) = R- *(~) (a), (57) 

where / and R are new constants found from continuity condition of the total wave function </? t otai and its derivative 
at the internal turning point at P ,i n t: 



we = Vtotal * (_) '' - ^otal * ( ~ 



(58) 



R = y total ^ ' - ytotal 



*(+) *("),' - <f(-) 

tt— dtp, int 



a— a tp , i n t 



On the basis of these solutions we obtain at the internal turning point a tp ,i n t the flux incident on the barrier, the 
flux reflected from it and the flux of mixing. The flux transmitted through the barrier was calculated at the external 
turning point a tP;ext . 



H. Penetrability and reflection: fully quantum approach versus semiclassical one 

Now we shall estimate through the method described above the coefficients of penetrability and reflection for the 
potential barrier with parameters A = 36, B = 12 A, A = 0.01 at different values of the energy of radiation E^. 
We shall compare the coefficient of penetrability obtained with the values given by the semiclassical method. In the 
semiclassical approach we shall consider two definitions of this coefficient: 

pWKB,(l) J_ pWKB,(2) 1 

Dcnctrabilitv rO. ' d 



penetrability q2 ' penetrability ^ ^ 2 ' (59) 

where 



(29 + 1/(29)^ 



(ext) 
tp 



6> = cxp J \V(a)-E\da. (60) 



(int) 
2 tp 
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TABLE I: The penetrability Pponctrabinty of the barrier and the duration r of the formation of the Universe defined by eq. (61) 
in the fully quantum and semiclassical approaches 



One can estimate also the duration of the formation of the Universe, using by definition (15) in Ref. [2]: 

r = 2 a tp . int p • (61) 

-^penetrability 

The results are presented in Tabl. I. In calculations the coefficients of penetrability, reflection and mixing are 
defined by eqs. (54), the fluxes by eqs. (52) (calculated P^^TMhty comc ide with -Pp^™ r abint y U P to tnc nrst ^ digits 
for energies in range < -E ra d < 2500). 
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From this table one can see that inside the entire range of energy, the fully quantum approach gives value for the 
coefficient of penetrability enough close to its value obtained by the semiclassical approach. This differs essentially 
from results in the non-stationary approach [2] . This difference could be explained by difference in a choice of the 
boundary condition, which is used in construction of the stationary solution of the wave function. 



I. The penetrability in the FRW-model with the Chaplygin gas 

In order to connect universe with dust and its accelerating stage, in Ref. [20] a new scenario with 
the Chaplygin gas was proposed. A quantum FRW-model with the Chaplygin gas has been constructed on the 
basis of equation of state instead of p (a) = p ra d(a)/3 (where p (a) is pressure) by the following (see also Refs. [5, 6]): 

Pch = -^-, (62) 

^Ch 

where A is positive constant and < a < 1. In particular, for the standard Chaplygin gas we have a = 1. Solution 
of equation of state (62) gives the following dependence of density on the scale factor: 

where B is a new constant of integration. Using parameter a, this model describes transition between stage, when 
Universe is filled with dust-like matter, and its accelerating expanding stage (through scenario of Chaplygin gas 
applied to cosmology, for details, see Refs. [7, 8, 20], also historical paper [13]). 

Let us combine expression for density which includes previous forms of matter and the Chaplygin gas in addition. 
At limit a — » eq. (63) transforms into the pdust component plus the p\ component. From such limit we find 

A = p A , B = p d u S t (64) 

and obtain the following generalized density: 

/ s l/(l+a) 

<•<«)- + im- < 65 » 

Now we have: 

o n ( / \ !/(!+«) 

2, ; o 71 " G ) 2 I „ , Pdust \ , Prad 



^'-t k + ^r°- (66) 



After quantization we obtain the Wheeler-De Witt equation 



^2 + ^ch (a) | <p(a) = £ rad <p{a), E r&d = (67) 



O \ 2 r, / N V(l + «) 

\ 7 2 ° 4 ( Fdust 



where 

For the Universe of closed type (at k = 1) at 87rG = M~ 2 = 1 we have (see eqs. (6)-(7) in Ref. [7]): 

/ Da <- \ !/(!+") 

U C h (a) = 36 a 2 - 12 a 4 (A + ^g^y) , ^rad = 12 /w (69) 

Let us expand the potential (69) close to arbitrary selected point a by powers of q = a — a and restrict ourselves to 
linear terms: 



^Ch (q) = V + V iq . 



(70) 
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FIG. 7: Cosmological potentials with and without Chaplygin gas: Left panel is for potential V{a) = 36 a 2 — 12 A a 4 with 
parameter A = 0.01 (turning point a tp = 17.320508 at zero energy -E ra( j = 0), Right panel is for potential (69) with parameters 
A = 0.01, pdust = 30, a — 0.5 (minimum of the hole is -93.579 and its coordinate is 1.6262, maximum of the barrier is 177.99 
and its coordinate is 5.6866). 



For coefficients Vq and V\ we find 

V = 
V x = 



Vq = Vch {a = a), 
dV Ch (a) 



da 



= 72 a + 12 a 3 <^ -4 A- 



{-«-^£t}-(a 



Pdust \- Q/(1+a) 



a 3(l+a) 



and eq. (67) has the form: 



After the change of variable 



cq. (72) becomes: 



After the new change 



we have 



From eqs. (73) and (75) we have: 



sp(q) + (Vq - E Iad + V x q) ip{q) = 0. 



^(0 



J ^rad — Vq V\ 

i i^i 2/3 "r^ 



C^(C)=o. 



-Erad — V) 

|Vi| 2 / 3 



v_ 

Wi\ 



E Ia d — Vq 

|Vil a /3 1^)2/3 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 



(77) 



Using such corrections after inclusion of the density component of the Chaplygin gas, we have calculated the wave 
function and on its basis the coefficients of penetrability, reflection and mixing by the formalism presented above. 
Now following the method of Sec. Ill A, we have defined the incident and reflected waves relatively to a new boundary 
which is located in the minimum of the hole in the internal region. Results are presented in Tabl. 3. One can see 
that penetrability changes up to 100 times, in such a coordinate, in dependence on the location of the boundary or 
in the internal turning point (for the same barrier shape and energy -E ra d)! This confirms that the coordinate where 
incident and reflected waves are defined has essential influence on estimation of the coefficients of penetrability and 
reflection. This result shows that the method proposed in the present paper has physical sense. In the next Tabl. 4, 
we demonstrate the fulfillment of the property (55) inside the entire energy range, which is calculated on the basis of 
the coefficients of penetrability, reflection and mixing obtained before. 
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V. MULTIPLE INTERNAL REFLECTIONS FULLY QUANTUM METHOD 



A. Passage to non-stationary WDW equation: motivations 



Tunneling is a pure quantum phenomenon characterized by the fact that a particle crosses through a classically- 
forbidden region of the barrier. By such a reason, the process of incidence of the particle on the barrier and its further 
tunneling and reflection are connected by unite cause-effect relation. So, the dynamical consideration of the tunneling 
process through cosmological barriers is a natural one. The rejection of the dynamical consideration of tunneling from 
quantum cosmology limits the possible connection between initial stage, when the wave is incident on 
the barrier, and next propagation of this wave. This leads to uncertainties in determination of penetrability 
and rates. According to quantum mechanics, a particle is a quantum object having properties both particle and wave. 
In the classically forbidden regions the wave properties of the studied object are evident. So, the wave description of 
tunneling is natural. 

So, we define a non-stationary generalization of WDW equation as 

j^-V«ff(a))#(a,T) = -i^¥(a,T), (78) 

where r is a new variable describing dynamics of evolution of the wave function being analog of time. According to 
quantum mechanics, the penetrability and reflection are stationary characteristics, and such characteristics, obtained 
in the following, are independent on the parameter r. Note that all these characteristics are solutions of stationary 
WDW equation, while non-stationary consideration of multiple packets moving along barrier gives clear understanding 
of the process. 

In order to give a basis to readers to estimate ability of the approach developed in this paper, let us consider results 
in [35] (see eq. (19)). Here was studied the non-stationary WDW equation 

i ^ - V eB (a)) *(a, r) = -i^ *(o, r) (79) 
with the potential for the closed FRW model with the included generalized Chaplygin gas. 



V eS (a)=3a*-^-\/A + ^ (80) 



After change of variable a new = a Q id \f\2 the non-stationary eq. (79) transforms into our eq. (78) simce the V e s 
potential is independent on the r variable (such a choice allows a correspondence between energy levels, convenient 
in comparative analysis). The potential (79) after such a transformation is shown in figs. 8. We shall analyze the 




Scale factor, a Scale factor, a 



FIG. 8: Behavior of the potential (80) after change a nC w = a id \/l2 at A = 0.001 and B — 0.001 (choice of parameters see in 
fig. 1, tables I and II in [35]): (a) shape of the barrier (V max = 223.52 at a = 42.322); (b) there is a little internal well close to 
zero (Vmin = -8.44 at a = 0.00581) 



behavior of the wave function. 



B. Tunneling of the packet through a barrier composed from arbitrary number of rectangular steps 



Now let us come to another more difficult problem, namely that a packet penetrating through the radial barrier 
of arbitrary shape in a cosmological problem. In order to apply the idea of multiple internal refections for study 
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the packet tunneling through the real barrier, we have to generalize the formalism of the multiple internal reflections 
presented above. We shall assume that the total potential has successfully been approximated by finite number N of 
rectangular steps: 



V(a) = { 



' V\, at a min < a < ai (region 1), 
V2, at ai < a < a 2 (region 2), 



(81) 



,V N , at ajv-i < a < a max (region N), 



where V. are constants (i = 1 . . . N). Let us assume that the packet starts to propagate outside inside the region with 
some arbitrary number M (for simplicity, we denote its left boundary <Im-i as a s tart) from the left of the barrier. We 
are interested in solutions for energies above that of the barrier while the solution for tunneling could be obtained 
after by change i£i — > ki. A general solution of the wave function (up to its normalization) has the following form: 



ip (a) = < 



-ik\a 



ai e"" 1 " + 0i e" 

at a m in < a < a\ (region 1), 

OLM-i e ikM ~ ia + (3 M -i e' lkM - ia 



at a.u-1 < a < «m-i (region M — 1), 



e ik M a + A R e 
at ttM-i < a < 



— ikM& 

a < «m (region M), 



(82) 



a M +i e lkM + ia + Pm+i e~ tkM + ia , 

at om < a < ajf+i (region M + 1), 



a n 



ik^-ia 



at ajv-2 < a < ajv-i (region N — 1), 



At e ikNa , at oat-i < a < a„ 



(region iV), 



where ctj and /3j are unknown amplitudes, At and are unknown amplitudes of transmission and reflection, 
ki = j-^2m(E — Vi) are complex wave numbers. We have fixed the normalization so that the modulus of the starting 
wave e lkMa equals to one. We look for a solution of such a problem by the approach of the multiple internal reflections. 

Let us consider the initial stage when the packet starts to propagate to the right in the region 
with number M. According to the method of the multiple internal reflections, propagation of the 
packet through the barrier is considered by steps of its propagation relatively to each boundary (see 
[10, 25, 27], for details). Each next step in such a consideration of propagation of the packet will 
be similar to the first 2N — 1 steps. From analysis of these steps recurrent relations are found for 



calculation of all unknown amplitudes Aj?\ 



A 



R 1 



and for arbitrary step n (for region with 



number j), summation of these amplitudes are calculated. We shall look for the unknown amplitudes, 
requiring the wave function and its derivative to be continuous at each boundary. We shall consider the coefficients 
T^, . . . and Rf , Rf ... as additional factors to amplitudes e ±lka . Here, the bottom index denotes the number of 
the region, upper (top) signs "+" and "— " denote directions of the wave to the right or to the left, correspondingly. 
To begin with, we calculate T^, T 2 ± . . . T^_ x and Rf, ■■■ Rn-i- 
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Analyzing all possible "paths" of the propagations of all possible packets inside the barrier and internal well, we 
obtain: 

E Anc = 1 + R M R M-1 + R M R M-1 ' R M R M-1 + — = 



n=l 



1 + E (#M #M-l) =" 

m=l 1 JX M JX A/I 



M ±L M-1 



E4 n) = EAt' 'fe 2 ^V 

71=1 v n=l ' v 

+ ' R N-1 R N-2 ' + — } = 

= (E^)-^-i' (84) 



E = + R M ' R M-1 + 



n=l 



+ ' Rm-1 ' R M-1 + ••• : 

/ +00 \ 



R M-1 ^Aj 



M „=l 



where 



1 _ fi j ^-1 



^7+1 = ^7+i + Tf+iRIT&i I 1 + E (^;+i^7) m ) = (85) 

x m=l ' 



RJ+1 ' £ 



+°° - n T+T+ 



/ l-OO _ \ 11 

= ^1(1+ e (^ + i^-r) - - y% 

v m=l ' 1 — Kj +1 Kj 



Choosing as starting points, the following: 

R-n-i = ^JV-H = ^Mi m = (86) 

we calculate the coefficients R%_ 2 ■ ■ ■ Rm-> r m+i ■ ■ ■ r n-i ano - ^m+i ■ ■ ■ ^jv-i- 

We shall consider propagation of all packets in the region with number M, to the left. Such packets are formed 
in result of all possible reflections from the right part of potential, starting from the boundary ajf. In the previous 
section to describe their reflection from the left boundary R to the right one, we used coefficient Rq . Now since we 
want to pass from simple boundary om-i to the left part of the potential well starting from this point up to a m ; n , we 
generalize the coefficient R^-i to R~m-i- The middle formula in (85) is applicable when we use eqs. (83) 
for definition of and Rf. Finally, we determine coefficients atj and j3j-. 



E«5 n) - + E(^7-i) ro ) = 

n—1 x m— 1 7 



1 — Rf Ri-i -^V~ 

(87) 

+°° / \ / +00 _ 

E = t+Ji + E (Rt R J-in R t = 

n—1 ^ m— 1 ' 



1-^7-1 ^ 



j 3 

3 
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the amplitudes of transmission and reflection: 

+00 , . +00 , . 

a t = E4 n) , a r = x;4"), 

+~ , . T+ +™o^ . . ( 88 ) 

n—l i- j n—1 

and coefficients T and i? describing penetration of the packet from the internal region outside and its reflection from 
the barrier 

rp _ fc W I A I 2 - I A I 2 T T - fcjV lf-+ I 2 

J-MIR—-, \ A T\ —Kline '^bar, J bar — ^ Jjy_i i _ . 

«M ' ' ' «M ' (89) 

RMIR = \A R \ = |^4inc I ' -Rbar, -Rbar = | -Rjtf | ■ 

Choosing a m ; n = 0, we assume full propagation of the packet through such a boundary (with no possible reflection) 
and we have Rq = — 1 (it could be interesting to analyze results with varying Rq). We use the test: 

\A T \ 2 + \A R \ 2 = 1 or Tmir + Rmir = 1. (90) 

km 

Now if energy of the packet is located below then height of one step with number m, then the following change 

k m ^i£, m (91) 

should be used for description of transition of this packet through such a barrier with its tunneling. In the case of a 
barrier consisting from two rectangular steps of arbitrary heights and widths we have already obtained coincidence 
between amplitudes calculated by method of MIR and the corresponding amplitudes found by standard approach of 
quantum mechanics up to first 15 digits. Even increasing the number of steps up to some thousands has the right 
accuracy to fulfill the property (90). 

In particular, we reconstruct completely the pictures of the probability and reflection presented in figs. 9 (a) and 
(b), figs. 10 (a) and (b), figs. 11 (b), but using such a standard technique. So, the result concerning the oscillating 
dependence of the penetrability on the position of the starting point a s t a rt in such figures is independent on the fully 
quantum method chosen for calculations. 

This is an important test which confirms reliability of the method MIR. So, we have obtained full coincidence 
between all amplitudes, calculated by method MIR and by standard approach of quantum mechanics. This is why 
we generalize the method MIR for description of tunneling of the packet through potential, consisting from arbitrary 
number of rectangular barriers and wells of arbitrary sizes. 



C. Results 



We have applied the above method to analyze the behavior of the packet tunneling through the barrier (80) (we 
used a n cw — >• \/l2a id). The first interesting result is a visible change of the penetrability on the displacement of the 
starting point a m ; n < a < a\, where we put the packet. Using the possibility of decreasing the width of intervals up 
to an enough small value (and choosing, for convenience, the width of each interval to be the same), we choose a m ; n 
as starting point (and denote it as a sta rt), from where the packet begins to propagate outside. We have analyzed 
how the position of such a point influences the penetrability. In fig. 9 (a) one can see that the penetrability strongly 
changes in dependence of a s tart for arbitrary values of energy of radiation -E rac ;: it has oscillating behavior. Difference 
between its minimums and maximums is minimal at a sta rt in the center of the well (i. e. its change tends to zero in 
the center of the well) , this difference increases with increasing value of a sta rt and achieves the maximum close to the 
turning point. With this result, we may conclude that exists a dependence of penetrability on the starting point a s tart 
of the packet. The coefficients of reflection, oscillations and penetration on the position of the starting point a s t a rt are 
presented in next figs. 9 (b), (c), (d) and have similar behavior. 

Usually, in cosmological quantum models the penetrability is determined by the barrier shape. In the non-stationary 
approach one can find papers where the role of the initial condition is analyzed in calculations of rates, penetrability 
etc. 3 But, the stationary limit does not give us any choice on which to work. We conclude: (a) the penetrability 



Such papers arc very rare and questions about dynamics have not been studied deeply. 
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FIG. 9: Dependencies of the coefficients of the penetrability Tb ar (a), reflection Rbnr (b), coefficient of oscillations K oac 
(c) and coefficient of penetration Tmir (d) in terms of the position of the starting point a sta rt for the energy E = 220 
(A = O.OOf , B = O.OOf, a max = 70. The total number of intervals is 2000, for all presented cases the achieved accuracy 
is |Tb ar + i?bar — 1| < 10~ 15 ). These figures clearly demonstrate oscillating (i.e. not constant) behavior of all considered 
coefficients on a s tart- 



should be connected with the initial condition (not only in non-stationary consideration, but also in the stationary 
one) , which determines position (coordinate) of maximum of the packet which begins to propagate outside (at initial 
time moment t = 0). (b) Even in the stationary consideration, the penetrability of the barrier should be determined 
in dependence on the initial condition. 

The first question is how much these results are reliable. In particular, how stable will such results be if we shift 
the external boundary outside? The results of such calculations are presented in fig. 10, where it is shown how the 
penetrability changes with a max (for clearness sake, we have fixed the starting point: a sta rt = 10). One can see that 
all calculations are well convergent, that confirms efficiency of the method of the multiple internal reflections. On the 
basis of such results we choose a max = 70 for further calculations. However, one can see that inclusion of the external 
region can change the coefficients of penetrability and penetration up to 2 times for the chosen energy level. 

The second question is how strong this affects the calculations of the penetrability. If it was small than, the 
semiclassical approaches would have enough good approximation. From figs. 9 it follows that the penetrability is 
not strongly changed in dependence on shift of the starting point. However, such small variations are connected 
with relatively small height of the barrier and depth of the well, while they would be not small at another choice of 
parameters (the coefficient of oscillation and penetration turn out to change at some definite energies of radiation, see 
below). So, this effect is supposed to be larger at increasing height of the barrier and depth of the well, and also for 
near-barrier energies (i. e. for energies comparable with the barrier height) and above-barrier energies of radiation). 

We have analyzed how these characteristics change in dependence on the energy of radiation. We did not expect 
the results that we got (see figs. 11). The coefficient of penetration has oscillations with peaks clearly 
shown. These peaks are separated by similar distances and could be considered as resonances in 
energy scale. So, by using the fully quantum approach we observed for the first time clear pictures of resonances 
which could be connected with some early unknown quasi-stationary states. At increasing energy of radiation the 
penetrability changes monotonously and determines a general tendency of change of the coefficient of penetration, 
while the coefficient of oscillations introduces the peaks. Now the reason of the presence of resonances has become 
clearer: oscillations of the packet inside the internal well produce them, while the possibility of the packet to penetrate 
through the barrier (described by the penetrability of the barrier) has no influence on them. In general, we observe 
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External boundary, a max External boundary, a m>x 

FIG. 10: Dependencies of the coefficients of penetrability (a), reflection (b), oscillations (c) and penetration (d) on the position 
of the external region, a max for the energy E = 223 (A — 0.001, B = 0.001). For all presented values we have achieved accuracy 
|T bar + i? bar -l| < l-HT 15 (the maximum number of intervals is 2000). 

134 resonant levels inside energy range S ra d = 0-200, and else 19 levels inside i? ra d = 200-223. 

In the last fig. 12 one can see that we have achieved |Tb ar + -Rbar — 1| < 10~ 15 inside whole region of changes of a s tart 
and a max (such data were used in the previous figs. 9 and 10). This is the accuracy of the method of the multiple 
internal reflections in obtaining Tb ar and i?bar- 

D. The fully quantum penetrability versus semiclassical one in cosmology: a quick comparison 

Does the penetrability, determined according to the semiclassical theory by a shape of the barrier between two 
turning points, give exhaustive answers and the best estimations of rates of evolution of universe? If we look at 
figs. 9 (a), we shall see that this is not the case. The penetrability is depended on the position (coordinate) of 
maximum of the packet which begins to propagate outside at time moment t — 0. So, the penetrability should be 
a function of some parameters of the packet at its start. For the first time, it has been demonstrated the difference 
between the fully quantum approach and the semiclassical However, let us perform a general analysis. 

(1) If we wanted to check the semiclassical approach, we should miss some of the parameters. One 
can use test of T + R = 1 (where T and R are the penetrability through the barrier and reflection from 
it). But, note that the semiclassical approximation neglects the reflected waves in quantum mechanics 
(see [21], eq. (46.10), p. 205, p. 221 222). Therefore, we cannot use the test above for checking T in 
the semiclassical theory. 

(2) If we would like to determine the reflection coefficient, then we should find a more accurate semiclassical 
approximation (in order to take into account both decreasing and increasing components of the wave function in 
the tunneling region). In such a case, we shall face another problem, namely the presence of a non-zero interference 
between the incident and reflected waves. Now the relation T + R = 1 cannot be used as test, and one needs to 
take the third component M of interference into account (see [31]). If we improperly separate the exactly known full 
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FIG. 11: Dependencies of the coefficient of the coefficient of penetration Tmir (a), the coefficient of the penetrability Tb ar (b), 
coefficient of oscillations K OS c (c) and difference E ICS , next — -Ercs, previous between two closest energy peaks (d) on the E la .d energy 
(we have choose: A = 0.001 and B = 0.001, a star t = 10, a max = 70, number of intervals inside the scale axis a 1000, number 
of intervals of energy 100000). Inside the energy region _E ra d = 200 — 223 we observe 19 resonant peaks in the dependencies of 
coefficients Tmir and K OEC while the penetrability increases monotonously with increasing the i? ra d energy. 
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FIG. 12: Accuracy of the obtained penetrability Tb ar and reflection 7?b ar for the energy E — 220 used in previous figs. 9 and 10. 
As a test, we calculate Tb ar + Rhai — 1 in dependence on the position of the starting point a star t (a) and the external boundary 
a max (b) (A = 0.001, B = 0.001, total number of intervals is 2000). 



wave function in the incident and reflected waves 4 , the interference component should increase without limit. In such 
a case, the penetrability and reflection can freely exceed unit and increase without limit. What is now the general 



4 However, the semiclassical approaches have no apparatus for such an analysis. 
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meaning of the penetrability? 

(3) We shall give only some examples from quantum mechanics, (i) If we consider two-dimensional penetration 
of the packet through the simplest rectangle barrier (with finite size), we shall see that the penetrability is directly 
dependent on direction of tunneling of the packet. So, the penetrability is not a single value but a function, (ii) 
If we consider one-dimensional tunneling of the packet through the simplest rectangular barrier, we shall obtain 
"interference picture" of its amplitude in the transmitted region, which is dependent on time and space coordinates 
and is an exact analytical solution. Of course, the stationary part of such a result exactly coincides with well known 
stationary solutions [27]. 

(4) A tunneling boundary condition [46, 47] seems to be natural and clear, where the wave function should represent 
an outgoing wave at large scale factor a. However, is such a wave free? In contrast to problems of quantum atomic 
and nuclear physics, in cosmology we deal with potentials, which modules increase with increasing the scale factor a 
(their gradients increase, which have sense of force acting on the wave). Therefore, in quantum cosmology we should 
define the boundary condition on the basis of the waves propagating inside strong fields (see [31]). 

These points destroy the semiclassical basis of the cosmological models. Now the statement concerning reliability 
of the semiclassical approach become a question of " faith" (note that this is widespread [31, 32]). The semiclassical 
approach could be compared with "black box", where deeper and more detailed information about the dynamics of 
the universe is hidden. 



VI. CONCLUSIONS AND PERSPECTIVES 

In this Chapter the closed Friedmann-Robertson-Walker model with quantization in the presence of a positive 
cosmological constant and radiation was studied. We have solved it numerically and have determined the tunneling 
probability for the birth of an asymptotically de-Sitter, inflationary Universe as a function of the radiation energy. 
Note the following. 

1 . A fully quantum definition of the wave which propagates inside strong field and which interact minimally with 
them, has been formulated for the first time, and approach for its determination has been constructed. 

2. A new stationary approach for the determination of the incident, reflected and transmitted waves relatively to 
the barrier has been constructed. The tunneling boundary condition has been corrected. 

3. A quantum stationary method of determination of coefficients of penetrability and reflection relatively to the 
barrier with analysis of uniqueness of solution has been developed, where for the first time non-zero interference 
between the incident and reflected waves has been taken into account and for its estimation the coefficient of 
mixing has been introduced. 

4. In this chapter the a development of the method of multiple internal reflections is presented (see Rcfs. [10, 25- 
28, 36], also Refs. [3, 16, 34]). When the barrier is composed from arbitrary number n of rectangular potential 
steps, the exact analytical solutions for amplitudes of the wave function, the penetrability Tb ar through the 
barrier and the reflection i?bar from it are found. At n — > oo these solutions can be considered as exact limits 
for potential with the barrier and well of arbitrary shapes. 

In such a quantum approach the penetrability of the barrier for the studied quantum cosmological model with 
parameters A = 36, B = 12 A (A = 0.01) has been estimated with a comparison with results of other known methods. 
Note the following. 

1. The coefficient of reflection from the barrier in the internal region has been determined. According to cal- 
culations, this coefficient is different visibly from unity at the energy of radiation close enough to the barrier 
height. 

2. The modulus of the coefficient of mixing is very small (it is less 10~ 19 ). This points out that there is no 
interference between the found incident and reflected waves close to the internal turning point. 

3. On the basis of the calculated coefficients we reconstruct a property (55) inside the whole studied range of 
energy of radiation (see Fig. 12). 

4. The probability of penetration of the packet from the internal well outside with its tunneling 
through the barrier of arbitrary shape is determined. We call such coefficient as coefficient 
of penetration. This coefficient is separated on the penetrability and a new coefficient, which 
characterizes oscillating behavior of the packet inside the internal well and is called coefficient 
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of oscillation. The formula found, seems to be the fully quantum analogue of the semiclassical 
formula of T width of decay in quasistationary state proposed in Ref. [18]. Here, the coefficient 
of oscillations is the fully quantum analogue for the semiclassical F factor of formation and the coefficient of 
penetration is analogue for the semiclassical T width. 

5. The penetrability of the barrier visibly changes in dependence of the position of the starting point i? s tart inside 
the internal well, where the packet begins to propagate (see figs. 9). We note the following peculiarities: the 
penetrability has oscillating behavior, difference between its minimums and maximums is minimal at i? s tart in 
the center of the well, with increasing i? star t this difference increases achieving to maximum near the turning 
point. The coefficients of reflection, oscillations and penetration have similar behavior. We achieve coincidence 
(up to the first 15 digits) between the amplitudes of the wave function obtained by such a method, and the 
corresponding amplitudes obtained by the standard approach of quantum mechanics (see Appendix B in [33] 
where solutions for amplitudes were calculated in general quantum decay problem). This confirms that this 
result does not depend on a choice of the fully quantum method applied for calculations. Such a peculiarity is 
shown in the fully quantum considerations and it is hidden after imposing the semiclassical restrictions. 

6. In the non-stationary and stationary considerations the penetrability of the barrier should be connected with 
the initial condition. We suggest that a possible introduction of the initial condition into the known stationary 
semiclassical models can change the obtained results. 

7. If one takes into account the external tail of the barrier, the penetrability is visibly changed. For example, the 
penetrability is changed up to 2 times (see figs. 10) for the barrier (8) with parameters A = 0.001 and B = 0.001 
(see figs. 8) at the energy of radiation £ ra( j = 223. If one increases the external boundary a max , all amplitudes 
and coefficients are convergent. This confirms efficiency of the developed method. 

8. The coefficient of penetration has oscillating dependence on the energy of radiation. Here, peaks are clearly 
shown. They are localized at similar distances (see figs. 11). So, for the first time we have obtained in the fully 
quantum approach a clear and stable picture of resonances, which indicate the presence of some early unknown 
quasistationary states. If the energy of radiation increases, the penetrability is monotonously changed. It 
describes a general tendency of behavior of the coefficient of penetration, while the coefficient of oscillations 
gives peaks. Now the reason of existence of resonances becomes clear: oscillations of the packet inside the 
internal well give rise to them. In particular, we establish 134 such resonant levels inside range i? ra( j = 0-223 
for the barrier (8) with parameters A = 0.001 and B = 0.001. 

9. A dependence of the penetrability on the starting point has maxima and minima. This allows to predict some 
definite initial values of the scale factor, when the universe begins to expand. Such initial data is direct result 
of quantization of the cosmological model. 

10. The modulus of the wave function in the internal and external regions has minima and maxima which were 
clearly established in [29, 31]. This indicates, in terms of values of the scale factor, where the probable 
"appearance" of the universe is the maximal or minimal. So, the radius of the universe during its ex- 
pansion changes not continuously, but consequently passes through definite discrete values connected with these 
maxima. It follows that space-time of universe on the first stage after quantization seems to be rather discrete 
than continuous. According to results [29, 31], difference between maxima and minima is slowly smoothed with 
increasing of the scale factor a. In this way, we obtain the continuous structure of the space-time at latter times. 
The discontinuity of space-time is direct result of quantization of cosmological model. This new phenomenon 
is the most strongly shown on the first stage of expansion and disappears after imposition of the semiclassical 
approximations. 
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